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We propose a novel class of temperature-tunable semiconductor metamaterials that exhibit nega¬ 
tive refraction in the terahertz spectral range. These metamaterials are based on doped semiconduc¬ 
tor superlattices with ultrathin barriers of about 1 nm thickness. Due to the tunnel transparency of 
the barriers, layers of the superlattice cannot be considered as isolated and, therefore, the classical 
homogenization approach is inapplicable. We develop a theory of quantum homogenization which 
is based on the Kubo formula for conductivity. The proposed approach takes into account the wave 
functions of the carriers, their distribution function and energy spectrum. We show that the com¬ 
ponents of the dielectric tensor of the semiconductor metamaterial can be efficiently manipulated 
by external temperature and a topological transition from the dielectric to hyperbolic regime of 
metamaterial can be observed at room temperature. Using a GaAs/Alo.sGao.yAs superlattice slab 
as an example, we provide a numerical simulation of an experiment which shows that the topological 
transition can be observed in the reflection spectrum from the slab. 


I. INTRODUCTION 

Hyperbolic metamaterials (HMMs) are one of the 
fastest developing branches of modern opticsii^— The di¬ 
electric function of HMMs is described by a tensor with 
two different components corresponding to the directions 
along (e||) and across (ex) the optical axis: 

(e± 0 0 \ 

e= 0 ex 0 . (1) 

Vo 0 eii/ 

Depending on the sign of these components, the crys¬ 
tal represents a dielectric medium (ex > 0, e|| > 0), 
a metal (ex < 0, ey < 0) or a hyperbolic metamaterial 
(exe|| < 0). For HMMs the shape of equal-frequency sur¬ 
face in k-space represents a one- or two-sheet hyperboloid 
depending on the signature of permittivity tensor.— This 
results in a singularity of the photon density of states and 
explains the unique optical properties of HMMs— 

In tunable metamaterials, the dielectric and magnetic 
responses and, therefore, the shape of equal-frequency 
surface can be manipulated by external influences such 
as, for example, DC magnetic field^Ti^ temperaturer^rJ^ 
femtosecond light pulses application of a volt- 

ageji^^— illumination)etc. In highly tunable meta¬ 
materials, the signature of the permittivity tensor can be 
switched, changing the equal-frequency surface topology. 
This phenomenon is called topological transitioni^ It has 
been experimentally observed is some systems— 

Fabrication of tunable THz metamaterials is an im¬ 
portant problem because of numerous potential appli¬ 
cations in far-held subwavelength imaging)^ enhanced 
nonlinearitieS)^ nanoscale wave guiding and strong light 
conhnementi^ Some realizations of tunable THz meta¬ 
materials were presented, for example, in Refs. [2^ and 
[^ . The majority of them represent an array of resonators 
whose capacity and/or inductivity is changed by external 
inhuencesi^ 


Here we propose a new concept of an ultra homo¬ 
geneous temperature tunable metamaterial based on a 
semiconductor superlattice for THz applications. Here, 
the term ultra homogeneous implies that the superlat¬ 
tice consists of coupled quantum wells separated by thin 
(~ I nm) tunnel-transparent barriers. Superlattices with 
barriers of such a thickness are widely used, for exam¬ 
ple, in quantum cascade lasersj^ In this case, in con¬ 
trast to a super lattice with thick barriers quantum 
effects are particularly relevant and, therefore, it is incor¬ 
rect to describe the dielectric function of each layer sep¬ 
arately and then apply the homogenization procedure. 
Therefore, another approach, which takes into account 
the wave functions of the carriers, their energy spectrum 
modified by the superlattice potential and the carrier dis¬ 
tribution function, should be used. We discuss the theory 
of proper approximation (quantum homogenization) fur¬ 
ther in Sec. H. 

It has been shown that highly doped semiconductor su¬ 
perlattice can exhibit properties of a hyperbolic medium 
at infrared frequencies The frequency range of the 
hyperbolic regime are defined by plasma frequency, which 
depends on the free carrier concentration. Free carrier 
concentration in semiconductors is extremely sensitive 
to the temperature, in contrast to dielectrics and met¬ 
als. For example, in the vicinity of the donor activa¬ 
tion temperature, it can change by several ordersWe 
have shown that high temperature sensitivity of plasma 
frequency in a semiconductor metamaterial can be ex¬ 
ploited for the efficient tuning of the metamaterial’s op¬ 
tical properties in the THz region and stimulation of the 
topological transition from the dielectric to hyperbolic 
regime. 

The paper is organized as follows. In Sec. H we de¬ 
velop a quantum homogenization theory and derive the 
main equations for the effective permittivity tensor. In 
Secs. HI and IV, we analyze the band structure and effec¬ 
tive dielectric function of a Te-doped GaAs/Alo.sGao.yAs 
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FIG. 1. (a) Conduction band profile of n-doped superlattice. Blue shaded areas represent the energy of the minibands. 
Green shaded area corresponds to the donor band. Thicknesses of the well and the barrier are denoted as di and d 2 , 
respectively, V is the height of the barrier, g{E) is a donor distribution function with standard deviation A and a maximum 
at Ed- (b) Electron energy dispersion in GaAs/Alo.sGao.rAs superlattice with following parameters: di = 10 mn, d 2 — I mn, 
V= 0.26 eV. 


superlattice depending on the temperature and frequency 
of the electromagnetic field. In Sec. V, numerical simula¬ 
tion of an experiment on the measurement of the reflec¬ 
tion spectrum from the superlattice slab is performed. 
In Sec. VI we discuss a figure of merit and tunability 
of semiconductor metamaterials. Finally, in Sec. VII we 
summarize our major results. 


II. MODEL 

A. Quantum homogenization 

Within the effective medium approximation, a multi¬ 
layered structure with layer permittivities £i and layer 
thicknesses di can be considered as uniaxial optical crys¬ 
tal with permittivity tensor © whose principle compo¬ 
nents are determined as 

— = -^—, e|| = d = ^di. ( 2 ) 

i i i 

As we have mentioned in the introduction, these formulas 
are inapplicable when the thickness of the layers is com¬ 
parable with electron wavelength and, therefore, quan¬ 
tum mechanics laws become relevant. We consider a 
more accurate approach based on the Kubo formulai^ 
It takes into account the distribution function of the car¬ 
riers, their wave functions and spectrum modified by the 
superlattice potential: 


£a{uj) = 


1 - 




0j{uj + fy) 


471"^ , , 

H-(3) 

UJ 


Here and in what follows, the index a = ||,_L corre¬ 
sponds to the directions along and across the optical axis 
of the metamaterial. Parameter is a permittivity of 
the lattice without free carriers, 7 is inverse momentum 
relaxation time of the carriers which is supposed to be 
isotropic for simplicity. 

The first term interprets classical Drude-Lorentz for¬ 
mula. One can see that implementation of a super¬ 
lattice in a semiconductor makes its plasma frequency 
anisotropic and we can distinguish plasma frequencies 
along (r 2 ||) and across (O^) the optical axis: 




47re^ 2 

e°“ (2Trh)^ 





dpi 


d^p. 


( 4 ) 


Here is the carrier energy in the i-th miniband which 
depends on the momentum p, f{E,pi^T) is the Fermi- 
Dirac distribution function, p is the chemical potential, 
T is the temperature. The sum is over all the minibands. 
Here we neglect hole contribution into the plasma fre¬ 
quency because we will consider n-doped semiconductor 
structures. Equation (|3|) is similar to the classical defini¬ 
tion of plasma frequency: 




47rne^ 


e^m 


( 5 ) 


Indeed, the difference between Eq. 0 and Eq. m is 
that in Eq. © we just average the inverse effective 
anisotropic mass 1 /m* = d^E/dp^ with distribution 
function f{E,p,,T). 

Second term in Eq. © describes interband transitions: 


(Jq, 


27 y fff 

(27r?i)3 ^ J JJ ujij - UJ + ij twjji 


( 6 ) 
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FIG. 2. Color online, (a) Frequency dependence of real part of dielectric function along (blue line) and across (black line) 
the optical axe. Temperature T = 300 K. (b) Temperature dependence of real part of dielectric function along (blue line) 
and across (black line) the optical axis. Permittivity of the media without free carriers £°° is shown by the dashed line. 
Frequency cu/27r = 3 THz. 


Here = (z, p| Ja |j, p) is the matrix element of the 
current operator between the unperturbed eigenstates 
of the superlattice and hwij = Ei — Ej. 

In order to calculate Hq and r2|| we need to determine 
the energy spectrum Eiijp) and the chemical potential fi. 


B. Energy spectrum of carriers 

Let us consider a periodic semiconductor superlattice 
with period d consisting of a quantum well with thick¬ 
ness di and a barrier with thickness d 2 and height V 
[Fig. dja)]. Effective masses in the well and barrier we 
put equal to mi and m 2 , respectively. 

The total energy of an electron in a superlattice can be 
represented as the sum of two terms which correspond to 
electron motion across and along the optical axis: 

E^{p) = ^+MPll) ( 7 ) 

Here mj_ is the effective mass of electrons in the direction 
across the optical axis. We take it equal to the effective 
electron mass in a bulk material. The second term of the 
right-hand side describes energy dispersion of electrons 
in z-th miniband. It can be found from the dispersion 
equation. 

cos(p||d/?i) = cos{pidi/h) cos{p 2 d 2 /h)— 


C. Chemical potential 

Efficient temperature manipulation of free carrier con¬ 
centration can be realized in narrow-gap or in doped 
semiconductor structures. Temperature tuning in meta¬ 
materials based on narrow-gap semiconductors was par¬ 
tially analyzed in Ref. [s^ Here we consider the case 
of doped semiconductor structures on the example of a 
superlattice with quantum wells uniformly doped with 
shallow donors. 

In highly doped structures, wave functions of neigh¬ 
bour donors can overlap. This results in a shift of donor 
levels and formation of a donor band. In the case of a 
considerable shift, the donor band can overlap with the 
conduction band. This phenomenon is called the Mott 
transition and will be discussed further. 

The chemical potential p can be calculated from the 
electroneutrality condition^^. It states that free carrier 
concentration is equal to the concentration of ionized 
donors: 


^ I fip,p,T) 


2d^p 

{2'KKf‘ 


nd 


g{E)dE 

2e(M-s)/(feT) I 


(9) 


Here Ud is the full donor concentration, g{E) is a donor 
distribution function which can be approximated by a 
Gaussian^ with standard deviation A and maximum at 

EdML 


^ ( d lk\ ■ I d lk\ f , P 2 TOi\ 

- -sin{pidi/h)sm{p2d2/h) -^- 

2 \P 2 m 1 pim2j 

( 8 ) 

where pi = y^2mi£’(p||), p 2 = y^ 2 m 2 (E(p||) - V). This 
dispersion equation can be obtained from the Schrbdinger 
equation using the Floquet’s theorem. 


III. BAND STRUCTURE OF THE 
SUPERLATTICE 

The model described above is applicable for superlat¬ 
tices of various compounds and designs. As an example, 
let us consider a superlattice that consists of GaAs quan¬ 
tum wells and Alo. 3 Gao. 7 As barriers. Thickness of the 
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FIG. 3. Color online. Temperature dependence of plasma 
frequency along the optical axis (blue line) and across the 
optical axis (black line). Insets show the shape of equal- 
frequency surfaces in the dielectric and hyperbolic regimes. 

quantum well di and the barrier c ?2 we put equal to 10 
and 1 nm, respectively. We consider the high frequency 
permittivity in Eq. (I3|) to be isotropic (ejj“ = e“) and 
put it equal to lli^ The calculated band structure and 
electron dispersion are shown in Fig. [TJb). We consider 
that quantum wells are uniformly doped with Te donors 
with concentration = 1 x 10^® cm“®. The energies 
of Te donors in bulk GaAs and in the superlattice are 
slightly different due to quantum confinement that arises 
from the superlattice potential. We neglect this differ¬ 
ence and put Ed = 0.03 eV as in a bulk materiali^ The 
standard deviation of the donor distribution function A 
for such a doping level is about several hundredths of an 
electron-volt. In the structure under consideration we 
put A = 0.03 eV, which is in accordance with Refs. [3^ 
and [ 3 ^ 

The width of the first miniband is about several hun¬ 
dredths of an electron-volt, which is comparable with the 
gap between the minibands. So, both weak-coupling and 
tight-binding approximations are poorly applicable for 
the structure under consideration. 

For the chosen parameters, the overlap of donor 
miniband and first miniband is approximately equal to 
0.01 eV. This means that even at zero temperature there 
are free carriers in the conduction band and their Fermi’s 
energy is about 0.01 eV. This corresponds to a concen¬ 
tration about 10^® cm“®. A simple estimation of plasma 
frequency using Eq. ([S]) yields II/(27r) ~ 1 THz. 


IV. EFFECTIVE DIELECTRIC FUNCTION OF 
THE SUPERLATTICE 

We calculate the frequency and temperature depen¬ 
dencies of the permittivity tensor components using the 
quantum homogenization approach [Eq. (|T])]. The fre¬ 
quency dependence of e± and e|| at room temperature 
is shown in Fig. UKa). The average energy between mini¬ 


TM polarization 



FIG. 4. Golor online. Frequency dependence of reflectance 
at three different temperatures. Shaded areas correspond to 
frequency bands of hyperbolic dispersion. The inset shows 
the relative orientation of the dielectric function (e^ and ey), 
the electric-field vector (transverse magnetic), the layered 
structure and a spherical prism with permittivity Sp = 11. 
Angle of incidence is /3 = 40°. 


bands and, therefore, frequencies of interband transitions 
are about 0.05 eV [Fig.[T](a)], which corresponds to a fre¬ 
quency of 10 THz. So, the contribution of interband tran¬ 
sition into the dielectric function [second term in Eq. ([3|)] 
can be neglected at frequencies of 1 THz without any 
considerable precision losses. Thus, quantum homoge¬ 
nization predicts that, beyond the interband transition, 
the tensor components of a superlattice with thin lay¬ 
ers can be described within the classical Drude-Lorentz 
formula with plasma frequency described by Eq. Q. 

This results qualitatively differ from effective param¬ 
eters obtained within the classical homogenization pro¬ 
cedure which predicts a resonance behaviour of ey at a 
nonzero frequency. There is no contradiction here. Clas¬ 
sical homogenization implies that all layers are isolated 
from each other and that the carriers do not move from 
one layer into a neighbouring one. Qualitatively it is 
equivalent to the restoring force that obstructs the car¬ 
rier transport. This force results in an appearance of 
the resonance in ey. Quantum homogenization implies 
that the barriers are tunnel transparent and charges can 
move freely throughout the whole volume of the sample. 
Therefore, the dielectric function £y is similar to that of 
a metal but takes into account interband transitions. 

One can see from Fig. ^a) that there are three fre¬ 
quency regions which correspond to different forms of 
the equal-frequency surfaces in k-space: (i) at frequen¬ 
cies a;/27r > 4.1 THz the material behaves like a di¬ 
electric; (ii) at frequencies 2.2 THz < w/27r < 4.1 THz 
the form of the equal-frequency surface is a hyperboloid 
and the material exhibits optical properties of HMM; (iii) 
at a;/27r < 2.2 THz electromagnetic waves decay expo¬ 
nentially into the medium, similar to the behaviour in a 
metal. 

The temperature dependence of e± and e|| at a fre- 
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quency of 3 THz is shown in Fig. [2jb). One can see that 
the material behaves as a dielectric, HMM or a metal 
depending on the temperature. Dielectric dispersion can 
be realized at the temperature of liquid nitrogen, the hy¬ 
perbolic regime is achieved at room temperature. Thus, 
a topological transition for THz radiation can be realized 
at temperatures reasonable for an experiment. 

The temperature dependence of the permittivity ten¬ 
sor components is weak at low temperatures. This is 
explained by the fact that the main part of donor levels 
are noticeably separated from the conduction band bot¬ 
tom, thus the activation of electrons on such levels occurs 
at temperature kT ~ Ed. At higher temperatures when 
all donors are ionized, the dielectric functions e± and e|| 
tend to constant values. In our case it occurs at tem¬ 
peratures much higher than the melting point of GaAs 
which is 1511 K.^ 

Figure [3] represents a topological phase state diagram. 
Solid lines show the temperature dependence of the longi¬ 
tudinal f2|| and transversal 12plasma frequencies. These 
lines divide the plane of the figure into three regions. It 
follows from Eq. m that every region corresponds to 
the one of the possible regimes: dielectric, metal or hy¬ 
perbolic. The shapes of equal-frequency surfaces corre¬ 
sponding to each regimes are shown in the insets of Fig. [S] 

One can see that plasma frequencies 12|| and 12are not 
equal to zero at low temperatures. This is explained by 
the overlapping of the donor and first minibands, the so- 
called Mott transition^. A part of the carriers does not 
freeze-out at low temperatures and makes a contribution 
into the plasma frequencies. 
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FIG. 5. Color online. Simulation of reflectance spectrum 
from 10 fim GaAs/Alo.sGao.rAs superlattice for different in¬ 
cident angles /3. Incident wave has TM-polarization. Param¬ 
eters of the superlattice is mentioned in Sec. imi Tempera¬ 
ture is taken equal to 300 K. Inverse momentum relaxation 
time is: (a) 7 = 1 x 10^^ s“^; (b) 7 = 1 x 10^® s“^. Black 
solid line corresponds to the Brewster’s angle which is de¬ 
termined by Eq. GQl). 


V. SIMULATION OF EXPERIMENT 

The efficiency of temperature tuning and the presence 
of a topological transition in a metamaterial based on a 
semiconductor superlattice can be confirmed experimen¬ 
tally by the measurement of reflection spectrum from the 
metamaterial. Here we provide a numerical simulation 
of a possible experiment.The scheme of the experiment 
is shown in the inset of Fig. SI A plane electromagnetic 
wave of TM-polarization is supposed to be incident at 
the angle /3 on the 10 /rm GaAs/Alo.sGao.r superlattice 
grown on an undoped GaAs substrate. The inverse mo¬ 
mentum relaxation time 7 in Eq. m we put equal to 
3 X 10^^ s“^. For the sake of simplicity, we neglect 
temperature dependence of 7 . In order to get efficient 
excitation of the optical states with high wave vectors, 
we consider that the wave incident on the sample passes 
through an undoped GaAs spherical prism with permit¬ 
tivity Ep = 11 . 

The frequency dependence of the reflectance for differ¬ 
ent temperatures and a fixed incident angle /? = 40° is 
shown in Fig. S] Frequency intervals of the hyperbolic 
regime are marked by shadows. For the temperatures 
300 K and 500 K there is a minimum of reflection in the 
hyperbolic regime. The nature of the minimum can be 


explained in terms of Brewster’s angle and will be dis¬ 
cussed further. For the temperature 77 K the minimum 
is smeared out because the frequency width of hyperbolic 
regime is comparable with Drude relaxation constant 7 . 

The frequency interval corresponding to the hyperbolic 
regime strongly depends on the temperature. So, at 77 K 
the metamaterial exhibits properties of HMM at frequen¬ 
cies 1.2 THz < a;/27r < 1.7 THz. At room temperature 
this frequency interval is 2.5 THz < a;/27r < 4.2 THz. 
Temperature drastically affects the reflection coefficient. 
For example, increasing of the temperature from 77 K to 
300 K at the frequency 2.5 THz results in tenfold increase 
of reflection coefficient from 0.03 to 0.3. 

Figure [5] shows the reflectance at room temperature 
versus incident angle and wave frequency. Typical values 
of inverse momentum relaxation time 7 for semiconduc¬ 
tors lie in the interval between 10 ^^ s“^ and 10 ^° s“^ 
(e.g., see Refs. HI and 0). Subfigures (a) and (b) corre¬ 
spond to these limit cases. One can see from Fig. EJa) 
that in the case of low losses, dielectric, metal, and hy¬ 
perbolic regimes have well-distinct frequency boundaries 
which smear out with increasing of the losses [Fig.EJb)]. 

In the dielectric regime (a;/27r > 4.0 THz) we have 
total reflection for glancing incident angles (/3 > 80°). 
Total reflection takes place because Ep > ey > ex > 
0. At incident angles less than 80° (/3 < 80°) we have 
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nearly total transmission and a minimum of the reflection 
coefficient in the vicinity of /3 = 50°. High transmission is 
explained by low optical contrast between the sample and 
prism in the dielectric regime. The minimum corresponds 
to the Brewster’s angle as- For an uniaxial crystal as 
can be determined as^ 


( 10 ) 

Frequency dependence of Brewster’s angle as is shown 
in Figs. EKa) and (b) by black solid line. 

In the metal regime < 2.2 THz) we have nearly 

total reflection for all angles of incidence except the min¬ 
imum in the vicinity of /3 = 85° at low frequencies. The 
minimum corresponds to the excitation of the surface 
plasmon polariton mode in the slab. 

In the hyperbolic regime (2.2 THz < a;/27r < 4.0 THz), 
there are regions of high and low reflection. The high re¬ 
flection region appears because the HMM under consid¬ 
eration does not support optical states with small lateral 
(lying in plane of the interface) components of the wave 
vector. Low reflectance in the region for /3 from 20° to 
50° can be explained by the Brewster’s angle when the 
TM polarized wave is not reflected from the sample. The 
minimum of reflection j3 ^ 60° is explained by the 
excitation of the waveguide mode in HMM slabj ^°d^ 


ub = Re <1^ arcsin 


/ex£|| 


— e 


£_Le|| - £, 



VI. FOM VS TUNABILITY 

The optimal operation frequency for semiconductor 
tunable metamaterials can be estimated from the condi¬ 
tion of high tunability and low losses. The main absorp¬ 
tion mechanisms in semiconductors are free carrier ab¬ 
sorption, fundamental absorption, and absorption due to 
the optical phononsi^ Fundamental absorption is dom¬ 
inant in the optical range. Absorption due to optical 
phonons plays an important role in polar semiconduc¬ 
tors. It has resonance behaviour and is significant around 
the wavelengths 30-35 ^m. Free carrier absorption in¬ 
creases with wavelength as square of the wavelength, so 
these losses make the main contribution in the THz re¬ 
gion. Rough estimation using Eq. ([ 3 ]) shows that figure 
of merit (FOM) rj for a;/27r=3 THz and 7 = I x 10^° s“^ 
is 


Re(e) w 
Im(£) 7 


( 11 ) 


Therefore, high FOM is achievable for uj ^ 

From the other hand, high tunability of optical prop¬ 
erties is reached near or below the plasma frequency [see 
Eq. ©]: 


00 ■ 


( 12 ) 


Typical effective electron mass m* in semiconductors 
varies from ~ O.OImg for narrow-band semiconductors 


(InAs, InSb) to ^ O.Ime (Si, GaAs), permittivity e°° ^ 
10 , electron concentration varies in a wide range from 
vanishingly small values to ~ I 0 ^° cm~° for highly doped 
semiconductorsSimple estimation using Eq. m yields 
the value of plasma frequency ^ I x 10^^ Hz 

for free carrier concentration I x I0^° cm“° and m* ^ 
O.Ime. It corresponds to the wavelength Ap ^ 3 /rm. 
Therefore, semiconductors can potentially be used for 
the fabrication of tunable metamaterials for near infrared 
applications. However it is very difficult to reach ef¬ 
fective tunability for highly doped semiconductors be¬ 
cause of the considerable screening and strong random 
potentialj °°i'^° The simplest way to avoid these problems 
is to decrease the doping level to lO^^-IO^® cm“°. Such 
concentrations correspond to plasma frequencies H/27r ^ 
0.3 - I X I0^° Hz, i.e. to the wavelength Ap ^ 30 - 
100 ^m. 

Therefore, effective tunability and low losses in semi¬ 
conductor metamaterials can be reached in the frequency 
range 7 ^ w < Hp which corresponds to THz frequency 
band. 


VII. CONCLUSION 

In this work we proposed a new concept of an ultra ho¬ 
mogeneous temperature-tunable metamaterial based on 
a doped semiconductor superlattice. We have shown that 
the classical homogenization procedure is inapplicable for 
the description of the metamaterial in terms of effective 
parameters because of the tunnel transparency of the 
barriers separating the quantum wells and that quan¬ 
tum homogenization should be used. We developed the 
theory of quantum homogenization applied to semicon¬ 
ductor nanostructured metamaterials. It is based on the 
Kubo formula for conductivity and takes into account 
wave functions of the carriers, their energy spectrum and 
distribution function. 

We have shown that the components of the dielectric 
tensor the semiconductor metamaterial can be efficiently 
manipulated by external temperature. Efficient temper¬ 
ature tunability is a distinctive feature of semiconductors 
which is explained by the high sensitivity of free carrier 
concentration to the temperature. On the example of a 
GaAs/Alo.sGao.rAs semiconductor superlattice with Te- 
doped quantum wells (n^ = I x 10^® cm“°) we have 
shown that the temperature of the topological transition 
from the dielectric to hyperbolic regime ~ 300 K for the 
frequency ~ 4 THz. Numerical simulation shows that 
the topological transition can be detected in an exper¬ 
iment on the measurement of reflection spectrum from 
the superlattice slab. 

A significant advantage of semiconductor metamateri¬ 
als is the possibility of their direct integration into op- 
tolectronic devices and optical integrated circuits. More¬ 
over, semiconductor materials combine two important 
features. On one hand, the energy spectrum of the car¬ 
riers in semiconductor nanostructures can be precisely 
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tailored with quantum engineering technologies. On the 
other hand, there are many methods of dynamic control 
of the electron distribution function in semiconductors, 
which are well-developed and widely applied in nano- 
and optoelectronics. These are, for example, electrical 
injection, optical pumping, thermal excitation, electron 
heating by electric field, etc. The advantages mentioned 
earlier and the rich functionality of semiconductor meta¬ 
materials allow to consider them as important element of 
future optoelectronics. 
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